Reading molecular messages from the intersections of high-order harmonic spectra at 

different orientation angles 
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We investigate the orientation dependence of high-order harmonic generation (HHG) from 
with different internuclear distances irradiated by intense laser fields both numerically and analyti- 
cally. The calculated molecular HHG spectra are found to be sensitive to molecular axis orientation 
relative to incident laser field polarization and internuclear separation. In particular, the spectra 
calculated for different orientation angles demonstrate a kind of intersection, which is identified 
as arising due to intramolecular two-center interference in the HHG. The striking "intersection" 
phenomenon can be used to probe the molecular instantaneous structure. 
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High-order harmonic generation (HHG) from aligned 
molecules in strong laser fields of femtosecond duration 
has proven to be a powerful tool for resolving and con- 
trolling the processes in an ultrafast time scale. For in- 
stance, recent experiments showed that calibrating the 
molecular recollision electronic wave packet by a refer- 
ence atom, the HHG spectra can be used to image molec- 
ular or bitalQ, and measuring the interference minima 
in the HHG spectra]!, H, 0, the HHG can also be 
used to probe molecular instantaneous structure @, 0|- 
Moreover, the investigations on the molecular orbital 
tomography 0, 0, [ll], Il2l and the effects of two-center 
interference in the HHG[I1, Q M, M, EE OS El are 
leading valuable insights into the mechanism of atomic 
and molecular HHG. 

Nevertheless, theoretical studies demonstrated that in 
some cases, the molecular properties already enter the 
recollision electronic wave packet. Accordingly, the spec- 
tral amplitude of the molecular recollision electronic wave 
packet is largely different from its reference atom in some 
energy regions [2Q] . This implies that accurately calibrat- 
ing the molecular recollision electronic wave packet can 
be difficult for some species of molecules. In addition, 
it was revealed that both two-center interference and 
the interference between different recombination electron 
trajectories]!]]] are responsible for the suppressed har- 
monics at certain orders. As a result, the minima in har- 
monic spectra may shift as the laser intensity changes [22l|. 

In the present paper, numerically investigating the ori- 
entation dependence of the HHG from H^ with different 
internuclear distances, we find that the molecular HHG 
spectra are sensitive to molecular axis orientation rela- 
tive to incident laser field polarization as well as to the 
internuclear separation. In particular, the spectra calcu- 
lated for different orientation angles demonstrate a kind 
of intersection, which is identified as arising due to in- 
tramolecular two-center interference in the HHG. Com- 



pared to the interference-related minima in the HHG 
spectra, the intersections of the harmonic spectra are 
easy to identify in practice. This finding is advocated 
for promising applications as a prospective tool to probe 
the molecular structure and dynamics. 

Let us consider the Hamiltonian of H^~ as (the atom 
units oth = e = m e = 1 are used throughout this paper) 



H(t)=p 2 /2 + V(r)-r--E(t), 



(1) 



where V{r) is the Coulomb potential and E(i) is the 
external electric field. In the 2D case, the Coulomb 
potential is written as V(x, y) = ~ z + 

1 _ V ' W; % /0.5+(a:+fl/2)2+j/2 

where Z is the effective charge, and 



% /0.5+(x-fl/2) 2 +y 2 

R is the internuclear separation, (for Z = 1 and R = 2 
a.u., the ground state energy for H^~ reproduced here is 
I p = 1.11 a.u.). In this paper, we assume that the molec- 
ular axis is coincident with the x-axis and the external 
field E(t) = £ , esin(woi) is linearly polarized with an ori- 
entation angle of to the molecular axis. E and uq are 
the amplitude and angular frequency of the laser field, 
e is the unit vector along the laser polarization. Our 
calculation will be considered for 780 nm trapezoidally 
shaped laser pulses with a total duration of 10 optical 
cycles and linear ramps of three optical cycles. Numeri- 
cally, the Schrodinger equation with the above Hamilto- 
nian H (t) is solved by the operator-splitting method with 
2048 time steps per optical cycle. The numerical conver- 
gency is checked using a finer grid. The coherent part 
of the harmonic spectrum is obtained from the Fourier 
transformed dipole acceleration expectation value, and 
only the harmonics polarized parallel to the incoming 
field are considered [2]. 

To analytically investigate the molecular HHG, we use 
the Lewenstein model|21| and consider the acceleration 
effect of the bound potential [l], The time-dependent 
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Figure 1: Harmonic spectra of 2D with Z — 1, R — 2 a.u. 
at different laser intensities and orientation angles 9, obtained 
by the exact numerical calculation. 

dipole moment is given bv[l5l ll6T| 

x(t) = ij^dr{^)d* rec (p st - A(t))e- iS &"*.*^ 

xE(i - r) • d ion (p st - A(i - r)) + c.c, (2) 

where 5(p st ,t,r) = / t *_ T dt" [(p st - A(t")) 2 /2 + ^] is the 
semiclassical action, A(t) = — f E(t')dt' is the vector po- 
tential of the external field, and p s t = J. dt"A(t")/T 
is the canonical momentum corresponding to the station- 
ary value. d ton (p) = (p|f|0) = (27r)~ 3 / 2 drexp(-p ■ 
r)r(r|0) is the bound-free dipole transition matrix ele- 
ment between the molecular ground state |0) and the 
continuum |p) in the ionization step, and d rec (p) = 
( Pk |f|0) = (2^)- 3 / 2 drexp(- Pk • r)r(r|0) is that 
in the recombination step which considers the effect of 
the electron acceleration in the vicinity of the parent 
ion before recombination, i.e., the effective momentum 
Pk = y/p 2 + 2J p p/ |p[ is adopted to describe the acceler- 
ation effect [l], y, B EDI • 

The wave function of the valence orbital of the H^ 
molecule with lsa g symmetry, investigated here, is ex- 
pressed in the LCAO-MO approximation 

^ lscJg (r) = N lsag (r + R/2) + la (r - R/2)] , (3) 

where Ni s is the normalization factor, <fii s is the atomic 
Is orbital in the configuration space, and R is the vector 
between the two atomic cores of the molecule. Then, the 
dipole transition moment for H^ can be written as 

d lsCT » = 2iN[ s<Tg [-cos(p ■ R/2)di.(p) 

+ sin(p-R/2)^ ls (p)R/2], (4) 
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Figure 2: Harmonic spectra of 2D H^ with Z = f , R = 2 a.u. 
at different laser intensities and orientation angles 8, obtained 
by the exact numerical calculation. 

where di s (p) is the atomic dipole moment from the Is 
orbital, and <^i s (p) is the Is wave function in the mo- 
mentum space. 

In Eq. 4, The factor cos(p ■ R/2) represents two-center 
interference [23], and the term, which is proportional to 
the internuclear distance R, leads to the breakdown of 
translation invariance. In our calculations, this term is 
omitted according to Ref . 0, HH, [2(| . Then, we obtain 

dZf g (p) - N Uag [-2i cos(p • R/2)d ls (p)]. (5) 

In Fig. 1, we plot the harmonic spectra of 2D 
with Z = 1,R = 2 a.u. at different laser intensities 
and orientation angles 9, obtained by the exact numeri- 
cal calculation. The comparison between the black and 
red curves in each subpannel of Fig. 1 shows that the 
harmonic spectra have a broad region within the HHG 
plateau with highly suppressed harmonic emission rate. 
The center of the suppressed region shifts to higher har- 
monic order as the orientation angle 9 increases. The 
broad suppressed regions arise from the two-center in- 
terference effect in the molecular HHG [2, 20]. But the 
locations of the minima in the broad suppressed regions 
are difficult to identify in some cases. For example, for 
9 = 50°, in Fig. 1(c), the minimum is clear at the 67th 
order. Fig. 1(b) is not as clear as Fig. 1(c). 

It has been revealed in Ref. [12] that the interference 
minima in the molecular HHG spectra may shift as the 
laser intensity changes. In particular, according to the 
simple point-emitter model proposed in RefQ, for Hj 
with R = 2 a.u. and 9 = 50°, the predicted interference 
minimum is at the 51st order. Compared to the observa- 
tion of the 67th order in Fig. 1(c), a large shift of about 
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Figure 3: Harmonic spectra of 2D HJ" with Z = 1, R = 2 
a.u. at different laser intensities and orientation angles 9, 
calculated using the Lewenstein model. 

16 orders occurs here. Based on the above analyses and 
discussions, we expect that probing molecular structure 
using the pronounced interference-related minimum can 
not be applicable in some cases. 

However, as we can see in each subpannel of Fig. 1, 
the two harmonic spectra at different orientation angles 
9 demonstrate an intersection in the plateau region. For 
example, in Fig. 1(c), the black curve for 9 — 40° is lower 
from the 19th to the 49th order, while the red curve for 
9 = 50° is lower from the 49th to the 93rd order. The 
striking intersection of the two curves is at the 49th or- 
der. As the laser intensity changes, the intersection of 
the two curves is almost invariable, as indicated by the 
vertical solid line. The calculated harmonic spectra at 
other orientation angles 9 also show the similar phenom- 
ena as those revealed in Fig. 1. But the intersections 
of the harmonic curves change as the orientation angles 
change. For example, as presented in Fig. 2, the inter- 
section of two harmonic curves at 9 = 30° and 9 = 50° 
is at the 39th order (indicated by the vertical solid line). 
This is different from that in Fig. 1. 

Next, we concentrate on the physical mechanism be- 
hind these phenomena. In Fig. 3, with the same parame- 
ters as in Fig. 1, we plot the harmonic spectra calculated 
using the Lewenstein model Eq. 2 and the modified tran- 
sition dipole Eq. 5. One can see that the primary charac- 
teristics of the harmonic spectra in Fig. 3 are analogous 
with those of the corresponding curves in Fig. 1. For 
instance, in each subpannel of Fig. 3, the intersection of 
two harmonic spectra in the plateau region is at the 49th 
order, as indicated by the vertical solid curves. This is in 
agreement with that in Fig. 1. The parallelism between 



Figure 4: Function curves of S(nioo) and harmonic spectra 
of 2D Hj with Z — 1, R = 1.7 a.u. at different orientation 
angles 6. (a): Function curves; (b), (c) and (d): harmonic 
spectra obtained by the exact numerical calculation. The 
laser intensity used here is I — 8.5 x 10 14 W/cm 2 . 

the corresponding curves in Fig. 1 and Fig. 3 shows that 
the modified model, i.e., Eq. 2 with Eq. 5, is applicable 
here for the description of the molecular HHG, especially 
for the angle dependence of the HHG. This applicability 
is also consolidated below. 

From the expressions of Eq. 2 and Eq. 5, we con- 
jecture that the effects of two-center interference are re- 
sponsible for the intersections of the harmonic spectra at 
different orientation angles. This conjecture is verified 
by the following analyses. According to Ref.P, Q, the 
harmonic intensity F(Q = hujq,9) along the laser polar- 
ization can be written as F(Q,8) cx J7 4 |a(il, 9)D(p, 9)\ 2 , 
where a(f2, 9) is the spectral amplitude of the molecular 
recollision electronic wave packet, D(p, 9) — e ■ d"*°^(p), 
and fl = E p + I p ^. E p is the electronic kinetic energy. 
The spectral amplitude a(fl, 9) is closely related to the 
ionization process[H [2l|. The transition dipole D(p, 9) 
corresponds to the recombination process. Both a(£l,9) 
and D{p, 9) are alignment dependent [20|. Here, we focus 
on D(p, 9). The angle 9 in the expression of D(p,9) is 
contained in the cos function in Eq. 5. We extract the 
cos function from D(p, 9), and write its power as 

S(nuJo) = cos(pR/2cos9)cos(pR/2cos9), (6) 

with the dispersion relation uloq — p 2 /2, that considers 
the electron acceleration before recombination[l]. n is 
the harmonic order. 

Based on the results shown in Fig. 3, we expect that 
Eq. 6 gives a good description of the orientation de- 
pendence of the HHG from H^~. In Fig. 4 and Fig. 5, 
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Figure 5: Function curves of S(nioo) and harmonic spectra 
of 2D with Z = 1, R — 2.5 a.u. at different orientation 
angles 9. (a): Function curves; (b), (c) and (d): harmonic 
spectra obtained by the exact numerical calculation. The 
laser intensity used here is I — 8.5 x 10 14 W/cm 2 . 

we show the comparisons between the function curves of 
Eq. 6 and the corresponding harmonic spectra at varied 
orientation angles 9 and internuclear distances R. These 
comparisons demonstrate the applicability of Eq. 6 in the 
prediction of the intersection of two harmonic spectra at 
different orientation angles 9. 

Specifically, for R = 1.7 a.u., in Fig. 4(a), the formula 
predicts that the intersections of two harmonic spectra 
are at the 45th order for 9 = 30° and 9 = 40°, the 53rd 
order for 9 = 30° and 9 = 50°, and the 61st order for 
9 = 40° and 9 = 50°. In the numerical cases in Fig. 
4(b)-(d), the corresponding intersections are at the 47th 
order, the 57th order and the 65th order, respectively. 
For R — 2.5 a.u., in Fig. 5(a), the formula predicts that 
those are at the 29th order for 9 = 40° and 9 = 50°, the 
35th order for 9 = 40° and 9 = 60°, and the 43rd order for 
9 — 50° and 9 = 60°. In Fig. 5(b)-(d), the corresponding 
intersections are at the 37th order, the 39th order and 
the 47th order, respectively. The large difference in the 
case of 9 = 40° and 9 = 50°, is due to the breakdown of 
the dispersion relation, used in our simulation, in the low 
energy region [21, |2{|. For R = 2 a.u., the formula also 
gives a good prediction. 

In addition, one can see from Fig. 4 and Fig. 5 that, 
the interference-related minima are not distinct in the 
harmonic spectra. It should be mentioned that our nu- 
merical method and obtained HHG spectra are compara- 
ble to the previous work of Ref. 0, H, S| • However, in con- 



trary to the claim of Ref. @, @, Q , we find that the broad 
suppressed regions in the harmonic spectra arise from the 
effects of two-center interference in the HHG, while the 
locations of the interference-related minima in the HHG 
spectra predicted in Ref. [1, d, 0] are difficult to identify, 
as shown above. Notice that in Ref. 0, 0], a spectrum- 
smoothing procedure is used to help the identification of 
the interference minima. This procedure is not adopted 
in our analysis, since we think it somehow ambiguous 
and we expect a "direct" comparison of the numerical 
observations to the experimental measurements. 

There are two factors those could influence the posi- 
tion of the interference minimum. First, two-center in- 
terference occurs not only in the recombination process 
of the HHG, but also in the ionization process of the 
HHG (into intermediate continuum states). The interfer- 
ence factor cos(P • R/2) is also included in the transition 
dipole di on (p) in Eq. 2 that corresponds to the ioniza- 
tion process. As a result, the spectral amplitude a(Q, 9) 
of the molecular recollision electronic wave packet, which 
is responsible for the fine structure of the molecular HHG 
spectrum, may be affected by the interference [20]. Sec- 
ondly, besides the ground state, the first excited state can 
also contribute to the harmonic emission in the broad 
suppressed region of the molecular HHG spectrum [30]|. 
Accordingly, the interference pattern may be modulated 
by the population of the first excited state. 

However, the parallelism between the predictions of 
Eq. 6 and the numerical results, as discussed above, re- 
veals that the intersections of the HHG spectra could be 
less influenced by the two factors. Particularly, compared 
to the interference-related minima, the intersections are 
easier to identify. These suggest that we can measure 
the molecular bond length through the measurement of 
harmonic spectra at different orientation angles 9, using 
the following equation 

S(n s uj , 0\) = S(n s 0J , 9 2 ), (7) 

where n s is the harmonic order that corresponds to the 
intersection of two harmonic spectra at the orientation 
angles 9\ and 9i. 

In conclusion, we have shown that due to the effects of 
two-center interference in the HHG, the harmonic spec- 
tra of Hrj~ at different orientation angles demonstrate the 
striking intersections in the high-frequency plateau re- 
gion. The phenomena discussed here are general. They 
are expected to appear in other species of molecules. Our 
results can be useful for promising applications allowing 
to use the "intersection" phenomenon as a prospective 
tool to probe the molecular structure and dynamics. 

This work is supported in part by Hong Kong Bap- 
tist University and the Hong Kong Research Grants 
Council, and NNSF(No.l0725521), 973 research program 
No.2006CB921400, 2007CB814800. 



■5 



Liu_Jie@iapcm.ac.cn 

J. Itatani, J. Levesque, D. Zeidler, Hiromichi Niikura, H. [15 
Pepin, J. C. Kieffer, P. B. Corkum, and D. M. Villeneuve, 
Nature 432, 867(2004). [16 
M. Lein, N. Hay, R. Velotta, J. P. Marangos, and P. L. 
Knight, Phys. Rev. Lett 88, 183903, (2002). [17 
M. Lein, N. Hay, R. Velotta, J. P. Marangos, and P. L. [18 
Knight, Phys. Rev. A. 66, 023805 (2002). 
M. Lein, P. P. Corso, J. P. Marangos, and P. L. Knight, [19 
Phys. Rev. A. 67, 023819 (2003). 
G. Lagmago Kamta and A. D. Bandrauk, Phys. Rev. A [20 
71, 053407 (2005). 

T. Kanai, S. Minemoto, and H. Sakai, Nature 435, [21 
470(2005). 

C. Vozzi, F. Calegari, E. Benedetti, J. -P. Caumes, G. 
Sansone, S. Stagira, M. Nisoli, R. Torres, E. Heesel, N. [22 
Kajumba, J. P. Marangos, C. Altucci, and R. Velotta, [23 
Phys. Rev. Lett. 95, 153902 (2005). 

J. Itatani, D. Zeidler, J. Levesque, Michael Spanner, D. [24 
M. Villeneuve, and P. B. Corkum, Phys. Rev. Lett 94, 
123902, (2005). [25 
J. Levesque, D. Zeidler, J. P. Marangos, P. B. Corkum, [26 
and D. M. Villeneuve, Phys. Rev. Lett. 98, 183903 
(2007). [27 
S. Patchkovskii, Z. Zhao, T. Brabec, and D. M. Vil- [28 
leneuve, Phys. Rev. Lett 97, 123004, (2006). 
R. Torres, N. Kajumba, Jonathan G. Underwood, J. S. 
Robinson, S. Baker, J.W. G. Tisch, R. de Nalda, W. A. [29 
Bryan, R. Velotta, C. Altucci, I. C. E. Turcu, and J. P. 
Marangos, Phys. Rev. Lett 98, 203007, (2007). 
Van-Hoang Le, Anh-Thu Le, Rui-Hua Xie, and CD. Lin, [30 
Phys. Rev. A. 76, 013414 (2007) 

Anh-Thu Le, X. M. Tong, and C. D. Lin, Phys. Rev. A. 
73 041402(R) (2006). 

C. Figueira de Morisson Faria, Phys. Rev. A. 76, 



043407(2007). 

T. Kanai, S. Minemoto, and H. Sakai, Phys. Rev. Lett. 
98, 053002 (2007). 

M. F. Ciappina, C. C. Chirila and M. Lein, Phys. Rev. 
A. 75, 043405 (2007). 

N. Wagner et ai, Phys. Rev. A. 76, 061403(R) (2007). 
M. F. Ciappina, A. Becker and A. Jaron-Becker, Phys. 
Rev. A. 76, 063406 (2007). 

V. I. Usachenko, P. E. Pyak, and S. I. Chu, Laser Phys. 
16, 1326 (2006). 

Y. Chen, Y. Li, S. Yang, and J. Liu, Phys. Rev. A. 77, 
031402(R) (2008). 

M. Lewenstein, Ph. Balcou, M. Yu. Ivanov, Anne 
L'Huillier, and P. B. Corkum, Phys. Rev. A 49, 
2117(1994). 

Y. J. Chen and J. Liu, Phys. Rev. A. 77, 013410 (2008). 
J. Muth-Bohm, A. Becker, and F. H. M. Faisal, Phys. 
Rev. Lett. 85, 2280(4) (2000) 

J. Chen and S. G. Chen, Phys. Rev. A. 75, 041402(R) 
(2007). 

D. B. Milosevic, Phys. Rev. A. 74, 063404 (2006). 

W. Becker, J. Chen, S. G. Chen, and D. B. Milosevic, 

Phys. Rev. A. 76, 033403 (2007). 

P. B. Corkum, Phys. Rev. Lett. 71, 1994(1993). 

Zhangjin Chen, Toru Morishita, Anh-Thu Le, M. Wick- 

enhauser, X. M. Tong, and C. D. Lin, Phys. Rev. A. 74 

053405 (2006). 

Xibin Zhou, Robynne Lock, Wen Li, Nick Wagner, Mar- 
garet M. Murnane, and Henry C. Kapteyn, Phys. Rev. 
Lett. 100, 073902 (2008). 

Y. J. Chen, J. Liu, and Bambi Hu, High-order Har- 
monic Generation in a Time-integrated Quantum Tran- 
sition Picture, unpublished. 



